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Abstract 

We consider helical configurations of a cholesteric liquid crystal (CLC) sandwiched between 
two substrates with homogeneous director orientation favored at both confining plates. We study 
the CLC twist wavenumber q characterizing the helical structures in relation to the free twisting 
number qq which determines the equilibrium value of CLC pitch, Pq = 2iT/qQ. We investigate the 
instability mechanism underlying transitions between helical structures with different spiral half- 
turn numbers. Stability analysis shows that for equal finite anchoring strengths this mechanism can 
be dominated by in-plane director fiuctuations. In this case the metastable helical configurations 
are separated by the energy barriers and the transitions can be described as the director slippage 
through these barriers. We extend our analysis to the case of an asymmetric CLC cell in which the 
anchoring strengths at the two substrates are different. The asymmetry introduces two qualitatively 
novel effects: (a) the intervals of twist wavenumbers representing locally stable configurations with 
adjacent helix half-turn numbers are now separated by the instability gaps; and (b) sufficiently 
large asymmetry, when the difference between azimuthal anchoring extrapolation lengths exceeds 
the thickness of the cell, will suppress the jump-like behaviour of the twist wavenumber. 
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I. INTRODUCTION 



In equilibrium cholesteric phase molecules of a liquid crystal (LC) align on average along 
a local unit director n(r) that rotates in a helical fashion about a uniform twist axis Q. 
This tendency of cholesteric liquid crystals (CLC) to form helical twisting patterns is caused 
by the presence of anisotropic molecules with no mirror plane — so-called chiral molecules 
(see 01 for a recent review). 

The phenomenology of CLCs can be explained in terms of the Frank free energy density 

/4n] = ^|iri(V ■ n)2 + [n ■ V X n + go]' 

+ /Cs [n X (V X n)]^ - ^^24 div (ndivn + n x (V x n))| , (1) 

where Ki, K2, and -ft'24 are the splay, twist, bend and saddle-splay Frank elastic constants. 
As an immediate consequence of the broken mirror symmetry, the expression for the bulk 
free energy (PJ contains a chiral term proportional to the equilibrium value of the CLC twist 
wavenumber, go- 

The parameter go, which will be referred to as the free twist wavenumber or as the free 
twisting number, gives the equilibrium helical pitch Pq = 27r/go. For the twist axis directed 
along the z-direction, the director field fi = (cosgo-z, singo-2, 0) then defines the equilibrium 
configuration in an unbounded CLC. Periodicity of the spiral is given by the half-pitch, 
Pq/2, because fi and — n are equivalent in liquid crystals. 

Typically, the pitch Pq can vary from hundreds of nanometers to many microns or more, 
depending on the system. The macroscopic chiral parameter, h = go-ft'2, (and thus the pitch) 
is determined by microscopic intermolecular torques 0| and depends on the molecular 
chirality of CLC consistuent mesogens. The microscopic calculations of the chiral parameter 
are complicated as it is necessary to go beyond the mean-field approach and to take into 
account biaxial correlations j3| . Despite recent progress jl[ 0] , this problem has not been 
resolved completely yet. 

In this paper we are primarily concerned with orientational structures in planar CLC 
cells bounded by two parallel substrates. Director configurations in such cells are strongly 
affected by the anchoring conditions at the boundary surfaces which break the translational 
symmetry along the twisting axis. So, in general, the helical form of the director field will 
be distorted. 

Nevertheless, when the anchoring conditions are planar and out-of-plane deviations of 
the director are suppressed, it might be expected that the configurations still have the form 
of the ideal helical structure. But, by contrast with the case of unbounded CLCs, the helix 
twist wavenumber g will now differ from go. 

It has long been known that a mismatch between the equilibrium pitch Pq and the twist 
imposed by the boundary conditions may produce two metastable twisting states that are 
degenerate in energy and can be switched either way by applying an electric field 0. This 
bistability underlines the mode of operation of bistable liquid crystal devices — the so-called 
bistable twisted nematics — that have been attracted considerable attention over past few 
decades [8, a,l3,El|. 

More generally the metastable twisting states in CLC cells appear as a result of interplay 
between the bulk and the surface contributions to the free energy giving rise to multiple 
local minima of the energy. The purpose of this paper is to explore the multiple minima 
and their consequences. 
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The free twisting number go cind the anchoring energy are among the factors that govern 
the properties of the multiple minima representing metastable states. Specifically, varying 
go will change the twist wavenumber of the twisting state, g. This may result in sharp 
transitions between different branches of metastable states. The dependence of the twist 
wavenumber g on the free twisting number go is then discontinuous. As far as we are aware, 
attention was first drawn to this phenomenon by Reshetnyak et. al 0|. 

These discontinuities are accompanied by a variety of physical manifestations which 
have been the subject of much recent important research. One such is a jump-like func- 
tional dependence of selective light transmission spectra on temperature as a result of a 
temperature-dependent cholesteric pitch, examined by Zink and Belyakov 0, 3|. More 
recently Belyakov et. al jisl EgI ] and Palto ^3 have discussed different mechanisms behind 
temperature variations of the pitch in CLC cells and hysteresis phenomena. 

In this paper we adapt a systematic approach and study the helical structures using 
stability analysis. This approach enables us to go beyond the previous work by relaxing 
a number of constraints. One of these requires anchoring to be sufficiently weak (where 
"sufficiently" will be discussed further below), so that the jumps may occur only due to 
transitions between the helical configurations which numbers of spiral half-turns differ by 
the unity jisl . 16l |. Noticeably, this assumption eliminates important class of the transitions 



that involve topologically equivalent structures with the half-turn numbers of the same 
parity. 

We shall also apply our theory to the case of non-identical confining plate and show that 
asymmetry in the anchoring properties of the bounding surfaces results in qualitatively new 
effects. Specifically, we find that sufficiently large asymmetry in anchoring strengths will 
suppress the jump-like behaviour of the twist wavenumber g when the free wavenumber go 
varies. 

The layout of the paper is as follows. General relations that determine the characteristics 
of the helical structures in CLC cells are given in Sec. |nj Then in Sec. IIIII we outline the 
procedures which we use to study stability of the director configurations. The stability 
analysis is performed for in-plane and out-of-plane fluctuations invariant with respect to in- 
plane translations. We study CLC cells with the strong anchoring conditions and the cases 
where at least one anchoring strength is finite. We formulate the stability conditions and 
the criterion for the stability of the helical structures to be solely governed by the in-plane 
director fluctuations. The expressions for the fluctuation static correlation functions are 
given. In Sec. IIVI we study the dependence of the twist wavenumber on the free twisting 
number. Finally, in Sec. |3 we present our results and make some concluding remarks. 
Details on some technical results are relegated to Appendix 1X1 



II. HELICAL STRUCTURES 
A. Energy 

We consider a CLC cell of thickness d sandwiched between two parallel plates that are 
normal to the z-axis: z = —d/2 and z = d/2. Anchoring conditions at both substrates 
are planar with the preferred orientation of CLC molecules at the lower and upper plates 
defined by the two vectors of easy orientation: e_ and e_|_. These vectors are given by 

e_ = e^;, e_|_ = cos A^e^: + sin A0ey, (2) 
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where A0 is the twist angle imposed by the boundary conditions. 

We shall also write the elastic free energy as a sum of the bulk and surface contributions: 

F[n]= [ Mn\dv+ V / iy,(n)ds (3) 

JV ^^^i Jz=vd/2 

and assume that both the polar and the azimuthal contributions to the anchoring energy 
H^(/(n) can be taken in the form of Rapini-Papoular potential [l8l |: 



wjn) = -4- 



z=ud/2 2 



l-(n,e.)^ l_(n,e,)^ , (4) 



z=ud/2 



where W^'^ and W^f^^ are the azimuthal and the polar anchoring strengths. 
The CLC helical director structures take the following spiral form 

no = cosm(2;) e^; + sinM(2) e^, u{z) = qz + (pQ, (5) 

where q is the twist (or pitch) wavenumber and 0o is the twist angle of the director in 
the middle of the cell. The configurations (jSj) can be obtained as a solution of the Euler- 
Lagrange equations for the free energy functional (0) provided the invariance with respect 
to translations in the x — y plane is unbroken. 

The translation invariant solutions can be complicated by the presence of the out-of- 
plane director deviations neglected in Eq. (0) and, in general, does not represent a helical 
structure. Using Eq. (0) is justified only for those configurations that are stable with respect 
to out-of-plane director fluctuations. The corresponding stability conditions will be derived 
in the next section. 

We can now substitute Eq. (0) into Eq. to obtain the following expression for the 
rescaled free energy per unit area of the director configuration (jSl): 

{2d/K2)F[no] = /[no] = {qd - q^df + 2 ^ w'^;^ sin^ u,, (6) 

v=±l 

where u^, is the angle between the vector of easy orientation and the director no at 
the plate z = I'd/2; and the dimensionless azimuthal anchoring energy parameter wlp^ is 
proportional to the ratio of the cell thickness, d, and the azimuthal extrapolation length. 



2K2 2&^' 



w\ ' = = — cosuy = (no,ei,) . (7) 

z=ud/2 



The energy © is of the well-known "smectic-like" form [3| and can be conveniently 
rewritten in terms of the following dimensionless parameters 

(3 = qd- Acf), po = qod - A0, a = 2(f)o - Acf) (8) 

by using the relations 

d{(l - Qo) = l3 - Po, 2ui, = vl3 + a. (9) 

Given the free twist parameter /3o the energy © is now a function of a and the twist 
parameter j3 which characterize the helical structure (0). The azimuthal angles of the director 
at the bounding surfaces, (py = u{yd/2), can be expressed in terms of the parameters (jH} 
and A0 as follows 

0+ = M+ + A0 = (/3 + a)/2 + A0, 0_ = u- = {-(3 + a)/2. (10) 
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B. Twist wavenumber and parity 



It is not difficult to show that in order for the configuration to be an extremal of the free 
energy these parameters need to satisfy the system of the following two equations: 

/3o = (3 + w'^hm{(3 + iya), u = ±1. (11) 

Equivalently, this system determines the extremals as stationary points of the energy (jH)) 
and can be derived from the condition that both energy derivatives with respect to a and 
f3 vanish. 

Eq. (|TT|l can now be used to relate the parameters a and /3 through the equation 

a = arctan[etan/?] + irk, e = —j—: j-r, (12) 

where k is the integer, A; G Z, that defines the parity of the configuration fi = (—1)'^. 
Indeed, substituting Eq. (fT2|l into Eq. (fTT| gives the relation between /?o and /? 

/5o = 7^(/?) =/5 + /ii(;J^+)sin(/3 + arctan[etan/?]), /i = (-1)'=, (13) 

that depends on k only through the parity. This remark also applies to the expression for 
the energy that after substituting the relation ()12j) into Eq. ® can be recast into the form 

fM = smv^Y + + - /i 5^ w[^^ cosv,, (14) 

u=±l 

where 

■^i^ = /5 + '^arctan[etan/9]. (15) 

In Sec. II VI we will find that there are different branches of metastable helical configurations. 
Each branch is characterized by the number of the spiral half-turns and fi is the parity of 
this number. For this reason, the integer k will be referred to as the half-turn number. 

Thus, we have classified the director structures by means of the parity /i and the dimen- 
sionless twist parameter /? that can be computed by solving the transcendental equation (fT^. 
Fig. HI illustrates the procedure of finding the roots of Eq. (fT^ in the /5 — 7 plane. 

In general, there are several roots represented by the intersection points of the horizontal 
line 7 = /3o and the curves 7 = 7±(/9). Each root corresponds to the director configuration 
which energy can be calculated from Eq. (|Tl|l . The equilibrium director structure is then 
determined by the solution of Eq. with the lowest energy. Other structures can be 
either metastable or unstable. 

C. Strong anchoring limit 

However, these results cannot be applied directly to the case of the strong anchoring 
limit, where W^'^^ 00 and the boundary condition requires the director at the substrate 
z = ud/2 to be parallel to the corresponding easy axis, n(z/c//2) || ey. 
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When the anchoring is strong at both substrates, it imposes the restriction on the values 
of q, so that q takes the values from a discrete set Q. This set represents the director 
configurations characterized by the parameter /3 and labelled by the half-turn number k 

(3 = qd-A(f) = TTk, keZ. (16) 

Substituting the values of (3 from Eq. (fTHIl into the first term on the right hand side of 
Eq. (jni) will define the equilibrium value of k as the integer that minimizes the distance 
between nk and /3o (= qod — Acf)). The step-like dependence of /3 on Pq for these equilibrium 
structures is depicted in Fig. |2fa). 

An experimentally important case concerns mixed boundary conditions in which the 
strong anchoring limit applies to the lower plate only, W^'^ oo. For brevity, this case will 
be referred to as the semi-strong anchoring. Now the relations |T2|l and |T3|l reduce to 

a = p + 2nk, (17) 

Po = P + wl^hm{2P) (18) 
and the energy of the helical structures ()14j) is now given by 

f{(3) = [w^^^ sm{2(3)] ' + 2^i+) sin^ (3. (19) 

Interestingly, in the semi-strong anchoring limit, the parity of half-turns, /i, does not enter 
either the energy p9|) or the relation ()18j) . 



III. STABILITY OF HELICAL STRUCTURES 

In this section we present the results on stability of the helical configurations (0). These 
results then will be used in Sec. II VI to eliminate unstable structures from consideration. We 
also give expressions for director correlation functions. These are in order to discuss the 
effects of director fluctuations. 

We begin with the general expression for the distorted director field 

n = cos 6 cos no + cos 6 sin ni + sin 9 n2, (nj, n^) = 6ij, (20) 

where the vectors ni and n2 are 

ni = — sin u{z) + cos u{z) ey, n2 = e^. (21) 

For small angles and 6 linearization of Eq. (j^Uj) gives the perturbed director field in the 
familiar form 

n^no + 5no, 5no = ni + 6* n2, -0=^^^, (22) 

where the angles and 6 describe in-plane and out-of-plane deviations of the director, 
respectively. 

Following standard procedure, we can now expand the free energy of the director field ()2U|) 
up to second order terms in the fluctuation field i/? and its derivatives 

F[n] ^ F[no] + F^'^i^], (23) 
F^'\'^]= [ fl,'\,p]dv+J2 I Wl'\^)As. (24) 

JV ^^^iJ z=ud/2 
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10 12 14 10 12 

P/7l P/7l 

FIG. 1: The curves representing the plot of the function 7+(/3) and 7_(/9) are shown as thick 
sohd and dashed lines, respectively. The points located at the intersection of the curves and the 
horizontal straight line J = Po give the roots of Eq. (|13|1 . The value of Pq is (10+ 1/2)tt (thin solid 
line) and is (10 + l/2)7r it w^f, (thin dotted lines). Two cases are illustrated: (a) w^^^ — '^^f,^ — 
= 7r/2; (b) wj,"^ = 10.0, = 2.19, = 7r/2 (see Eq. 



The second order variation of the free energy is a bihnear functional which represents 

the energy of the director fluctuations written in the harmonic (Gaussian) approximation. 
^From Eqs. (0)-® we obtain expressions for the densities that enter the fluctuation en- 
ergy 



2/f = K, ( Vi</. + V2^) ' + K2 {V,9 - V2 
+ [(Vo0)' + (Vo^)'] + qK, [29 Vo0 + qO'] 

- 2uK2i e Vi( 



>2 cos2u^ + w!:''^e'^ 



=i/d/2 



(25) 



(26) 



where Vj = (rij, V). 

In what follows we shall restrict our consideration to the case of fluctuations invariant with 
respect to in-plane translations, so that ip = ip{z). This assumption, although restricting 
applicability of our results, allows us to avoid complications introduced by inhomogeneity 
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FIG. 2: Dependence of the twist parameter (3 {= qd — Ac/)) on the free twisting parameter [5q 
(= q^d — Acf)) at wjp ^ = wjp'^^ = for various values of the dimensionless azimuthal anchoring 
parameter (= Wi^d/2K2): (a) strong anchoring hmit, W(f) — > oo, discussed in Sec. Ill (]| (b) 
= vr; (c) = tt/2; (d) = O.OSvr. Discussion of the cases (b)-(d) can be found in Sec. IIV Al 



of the helical structure In this case the fluctuation energy per unit area is 

rd/2 

z=ud/2 



l-d/2 



(27) 



where 5* is the area of the substrates. The operator K is the differential matrix operator 
that enters the linearized Euler-Lagrange equations for the director distribution pn|l . i.e. 



Kil) = 0. 



(28) 



The eigenvalues of K form the fluctuation spectrum pj| . The eigenvalues A can be computed 
together with the eigenmodes ij^^ by solving the boundary-value problem: 



z=ud/2 



A5 



0. 



(29) 
(30) 
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The expressions for K and Q'''^^ are given by 



K^d. 
K^d, 



W^-!^ cos2n. 



+ 1*0 '''' 



and Kq is the effective elastic constant 



K^ = K^-1K2{\~q^lq). (33) 

^From Eqs. ()31 |) -()32 |) the operators K and Q^^^ are both diagonal, so that the in-plane and 
out-of-plane fluctuations are statistically independent and can be treated separately. 



A. In-plane fluctuations 

1. Strong and semi-strong anchoring 

We begin with the limiting cases discussed in Sec. Ill CI and where at least one of the 
bounding surfaces imposes the strong anchoring boundary condition. This is the case when 
we have to use stability criterion related to the fluctuation spectrum which requires all the 
eigenvalues to be positive, A > 0, as to ensure the positive definiteness of the fluctuation 
energy 

It is not difficult to see that, for the strong azimuthal anchoring present at both substrates, 
the lowest eigenvalue is 

= K,{7r/df (34) 

and all the structures with the twist parameter (fTBj) are locally stable with respect to in-plane 
fluctuations. 

For the semi-strong anchoring with "* oo, this is no longer the case. The stability 
condition is now given by 

1 + 2wJ+^ cos(2/?) > 0. (35) 

/^From Eq. p8j) this condition requires the free twist parameter (5q to be an increasing 
function of the twist parameter p. It is derived in Appendix |X1 [see Eq. (jAllj) with w+ 
replaced by w^^'* cos(2/3)]. 



2. Weak anchoring 

We have a somewhat different situation when the azimuthal anchoring strength is finite 
at both substrates. In this case the stability conditions can be derived using an alternative 
procedure [l^ . The procedure involves two steps: (a) solving the linearized Euler-Lagrange 
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equations (j2Hl); and (b) substituting the general solution into the expression for the fluctu- 
ation energy (j2Zj). The last step gives the energy (j^Tj) expressed in terms of the integration 
constants, so that the stability conditions can be derived as conditions for this expression 
to be positive definite. 

Following this procedure, we can obtain the stability conditions for the helical structures 
characterized by the parity fi and the twist parameter jS related to the free twist parameter 
f3 through the relation The final result is 

Hf, = Af, + 2wJ+^wJ,"^ COSV+ cosv_ > 0, (36) 

An = ^i{w^^^ cosf+ + w^"'' cosf„) > 0, (37) 

where v± are defined in Eq. ()15|1 . These inequalities also follow immediately from the 
stabihty conditions ()A10|1 obtained in Appendix |X1 by putting w± = fJ^w^^^ cosv±. 

Violating either of Eqs. (pK|l - (jH7jl will result in instability caused by slippage of the director 
in the plane of the spiral. Such an instability cannot occur when the azimuthal anchoring 
is strong at both substrates. 



B. Out-of-plane fluctuations 



We now study stability of the helical structures with respect to the out-of-plane fluctua- 
tions. To this end we replace A with Ki{2/d)^X and rewrite the eigenvalue problem (j^ - lp]! 
for 6 in the following form: 



0, 



=±1 



r, = {qdfKjK^ = (/? + M){r^{l3 + A0) + 2r2(/3o - /5)), 



w, 



2Ki 



2L 



(38) 
(39) 
(40) 

(41) 



where t = 2z/d, Vi = Ki/ Ki and L^g'^ is the polar anchoring extrapolation length. 
The stability condition A > can now be readily written as follows 

4A^ + > 0, 

where A^ is the lowest eigenvalue of the problem (jnHj) -(|H ^ computed at = 0. 



(42) 



When the polar anchoring is strong at both substrates, 
known [see remark at the end of Appendix IX]: 



(±) 



A. 



7rV4. 



oo, the eigenvalue Am is 



(43) 



Otherwise, Km is below it/ 2 and can be computed as the root of the transcendental equation 
deduced in Appendix 1X1 [see Eq. (IA12|) ] 



- ft^) Sm2Km + Km{wl^^ + Wg ') COS 2/?^ = 



(44) 



where < < vr/2. 
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1. Strong anchoring 



In the strong anchoring hmit, Eq. p6|) imphes that the values of the twist parameter (3 
are quantized and do not depend on the free twist parameter [3q. Unstable configurations 
are characterized by twist wavenumbers violating the stability condition ()42|1 with Am given 
in Eq. ()4H|1 . These wavenumbers are described by the inequalities 

1/^- < 1/^ < /3 = qd/ir, ^0 = qod/TT, 

l/P± = -r2/3o ± [Mf + 2r2 - rg] . (45) 

These inequalities yield two different sets of unstable structures depending on the sign of 
the difference {2K2 — K3). For qo > these sets are given by 

2K2 > K3 : 

q>{n/d)\P+\ or q < -{n/d)0.\, (46) 
Ks > 2K2 : 

- (7r/rf)|;3+| < q < -(7r/c/)|/3_| at go > ti / {dr^Wr^ - 2r2. (47) 

Eq. (j^6|l shows that, when the energy cost of bend is relatively small, there are an infinite 
number of unstable configurations and the configuration loses its stability as the distance 
between its wavenumber q and go becomes sufficiently large. 

Otherwise, unstable configurations may appear only if the free wavenumber go exceeds 
its critical value given in Eq. ()47|1 . In this case the number of the unstable configurations 
is finite. From Eq. (jTfj) there is no unstable configurations for nematic liquid crystals with 
go = 0. This result has been previously reported in Ref. [2^ . 



2. Weak anchoring 

We now pass on to the case where the strengths of anchoring PF^^"* are not infinitely 
large. By contrast to the case of strong anchoring, the twist parameters (3 and /3o are 
now not independent. Rather we have the stationarity condition (|T3|l relating (3 and /Jo- 
in addition, if the polar anchoring is also not infinitely strong, the eigenvalue Am can be 
considerably reduced. 

In these circumstances, it is reasonable to approximate the left hand side of the stability 
condition ()42|) by its lower bound derived in the limit of weak polar anchoring, W^'^'^ 0, 
where Am vanish. Technically, the resulting condition 

r,(/?, ^l) = {(3 + A0) ( 2 r2 /?o(/?, ^A + (^3 - 2r2)(3 + r,^<P ) > (48) 

is sufficient but not necessary for stability. Thus, when the inequality ()48|) is satisfied, the 
structure will certainly be locally stable with respect to out-of-plane fluctuations whatever 
the polar anchoring is. 

Eq. (PH|) can now be used to study out-of-plane fluctuation induced instability of the 
helical structures which are otherwise stable with respect to in-plane fluctuations and thus 
meet the stabihty conditions (jSEI)-(ISZ|)- For Acj) = and positive twist wavenumbers, such 
instabilty may occur only if the doubled twist elastic constant exceeds the bend elastic 
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constant, 2K2 > K3, and the azimuthal anchoring energy is sufficiently large. In this case, 
however, can be made non-negative by increasing the value of the free twist parameter 
/3o- In other words, if the ratio of the cell thickness and the equilibrium CLC pitch is large 
enough to meet the condition PH|) we can neglect out-of-plane deviations of the director and 
use the "smectic-like" free energy 



C. Correlation functions 



Our calculations of the director fluctuation static correlation function {i/)[z)ij^{z')) use 
the relation [l^ 

{.Piz)i^iz')) = ^Giz,z'), (49) 



S 



where ks is the Boltzmann constant and T is the temperature. The Green function G{z, z') 
can be computed as the inverse of the operator K defined in Eq. by solving the 
boundary-value problem 



kG{z,z') =S{z-z')i, 



Q^''^G{z,z') 



z=ud/2 



0, 



(50) 
(51) 



where / is the identity matrix and the operator Q''"^ is given in Eq. (jH^ . Since the matrix 
operators K and Q'^'^^ are both diagonal, the correlation function is also diagonal: 



[tPizmzJj-l Q {e{z)9iz'))J- 



(52) 



Solving the boundary-value problem ()5m) - ()5H) yields the following expressions for the 
in-plane and the out-of-plane components of the correlation function (j52|) 



ksTd 



Wg ^ smK{l + 2z^/d) 



2SKikD{k 

+ /tcos/t(l + 2z^/d) w^g^^ sinK(l — 2zy/d) + kcosk(1 — 2zy/d) 



(53) 



(0(^)0(/)) 



ksTd 
2SK2H^ 



1 + /i(l + 2z^/d)w[ ■* cosf 



1 + /i(l — 2zy / d)w^^ cosf+ 



(54) 



where = —rq/4, z^ = min{z, z'}, z^ = max{z, z'}, D{k) and are defined in Eq. ()44|l 
and Eq. ()36|). respectively. The correlation functions diverge on approaching the boundary 
of the stability region. For out-of-plane fluctuations, the denominator of the expression 
for {6{z)6{z')) vanish in the limit of marginal stability where — Tq/4 = i^m and 

-D(k) — > D{nrn) = 0. Similarly, Eq. (jHUj) shows that, in the marginal stability limit for 



in-plane fluctuations, 
divergent. 



if^ goes to zero thus rendering the correlation function {(f){z)(f){z')) 
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FIG. 3: The free energy /+ (solid line) and /_ (dashed line) of stable configurations with the half- 
turn number, A;, between 10 and 14 as a function of /3o computed from Eq. (|14|) by using Eq. (|13|1 
for two values of the anchoring energy parameters: (a) Wff^ = tt/2; and (b) = tt. The intersection 
point of the branches with A; = 10 and /c = 11 at /3o = lO.Svr is indicated by circle. Squares mark 
energy of the structures at the transition point with /?o = llvr where the configuration with k = 10 
loses its stability. It is shown that the energy of the helical structure with A; = 10 (A: = 11) increases 
(decreases) as /?o varies from lO.Svr to llvr. 

IV. TRANSITIONS INDUCED BY FREE WAVENUMBER VARIATIONS 

In the previous section we have studied the stability of the CLC helical structures (0) with 
respect to both in-plane and out-of-plane fluctuations. We have found that the anchoring 
conditions play a crucial role in the calculations. In particular, cells with strong anchoring 
and those with what we have called semi-strong anchoring exhibit significantly different 
properties. 

In this section we concentrate on the weak anchoring cases. We have shown that in this 
case helical structures are characterized by the twist parameter, jS, and the half-turn parity, 
fi. These quantities are related to the free twist parameter, /?o, through the stationary point 
equation (fT^ . The structure responds to variations of the free wavenumber (and thus the 
free twist parameter) by changing its twist parameter. 

This change may render the initially equilibrium structure either metastable or unstable. 
When the anchoring is not infinitely strong and the free twist parameter is large enough 
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p/Ti; p/Ti; 

FIG. 4: The free energy /+ (solid line) and /_ (dashed line) of the configurations with the half-turn 
number, A;, between 10 and 14 as a function of (3 computed from Eq. (|14j) . Thin lines represent 
the energy of unstable configurations. Two cases are shown: (a) ^ = w^^"^ = = 7r/2; and 

(b) u;^ ^ = 10.0, w^^^ = 2.19 {a = 5.6 and = vr/2). Circles and squares label energy of the 
structures at [3q = lO.Svr and /3o = Utt shown in Fig. |2l^a). At /3o = lO.Svr, the structures with 
A; = 10 and A; = 11 are shown to be degenerate in energy and separated by the energy barrier A/. 

to meet the stability condition ()48|) . this instability is solely governed by in-plane director 
fluctuations and defines the mechanism dominating transformations of the director field. 
This mechanism is suppressed in the strong anchoring regime, where the structural transi- 
tions involve tilted configurations 0|, and can be described as director slippage through 
the energy barriers formed by the surface potentials. 

In this section our task is to study helical structure transformations as a function of 
the free twist wavenumber go for different anchoring conditions. Equivalently, we focus our 
attention on the dependence of (3 on (3q\ this can be thought of as a sort of dispersion relation. 
To this end we examine in more detail the consequences of the analytical results obtained 
in the previous sections. Sec. HT] and Sec. IIIII 
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FIG. 5: Dependence of (3 on /3o calculated at W- = 10.0 for various values of the parameters a 
(= 2w'f^wlj^^ /{w'f^ - wl^^)) and (see Eq. (a) = n {a = 11.8 and wl^^ = 3.705); (b) 

= tt/2 {a = 5.6 and wl^^^ = 2.19); (c) a = 1.0 (u^i^^ = w^c = 0.476); (d) a = 0.5 (u-i^^ = 0.25). 



A. Symmetric cells 

When the anchoring strengths at both substrates are equal, VF^""* = W^J"^"* = W^, the 
right hand side of Eq. (fT^ is /5 ± w^sinP {Wfj, = W(j,d/2K2) and v± = [3. In this case 
stability of the configurations is governed by Eq. (jHTjl which reduces to the simple inequality 
/X cos /3 > 0. 

It immediately follows that the values of (3 representing the locally stable structures of 
the parity /i ranged between {k — l/2)7r and {k + l/2)7r, where k is the even (odd) integer 
at /i = +1 (/i = —1). The integer k will be referred to as the half-turn number. The parity 
/i introduced in Sec. |n]is now shown to be the parity of the half-turn number: = (— l)'^. 

The intervals of (3 for the stable configurations [(fc — l/2)7r, (fc + l/2)7r] are now labelled by 
the half-turn number k. Since the function 7^ (/x = (— l)'^) monotonically increases on the 
interval characterized by the half-turn number fc, the value of /3o runs from {k — l/2)7r — 
to {k + l/2)7r + on this interval. As a result, for each half-turn number k, there is the 
monotonically increasing branch of the (3{(3q) curve. 

The branches with k ranged from 10 to 13 for different values of the dimensionless anchor- 
ing energy parameter are depicted in Figs. |2tb)-(d). We see that the /5o- dependence of 
/? will always be discontinuous provided the anchoring energy is not equal to zero. Fig. Efd) 
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shows that the jumps tend to disappear in the hmit of weak anchoring, where the azimuthal 
anchoring energy approaches zero, ^ 0. 

As we pointed out in Sec. Ill (^I for the case of strong anchoring, there are two equihbrium 
structures of the same energy at = (1/2 + m)7r. In Fig. Efa) the arrows indicate that the 
half-turn number of the equihbrium structure changes at these points. 

Similarly, when the anchoring is weak, < vr, and Po = (1/2 + m)7r, Eq. (fT!^ possess 
two different roots with k = m and k = m + 1 which are equally distant from /?o and are of 
equal energy. In Fig. Efa), the free energy (fT^ is shown as a function of Pq. It can be seen 
that the intersection points of the curves for different parities, (solid and dashed lines in 
Fig. El) are indeed at /3o = (1/2 + m)7i. The parity of the equilibrium configuration reverses 
as Po goes through the values (1/2 + m)7r. Fig.E^b) illustrates that this is also the case even 

if > 71. 

For m = 10, Fig. Efa) and Fig. ^a) show the initially equilibrium structure with the 
half-turn number k = 10 (solid line) becomes metastable as Po passes through the critical 
point Po = (1/2-1- 171)71 at which the structures with k = 10 and k = 11 (dashed line) are 
degenerate in energy. In Fig. Efa) and Fig. the structures at this point are indicated 
by circles. 

As is seen from the figures, relaxation to the new equilibrium state will require the jump- 
like change of the twist parameter p. In addition, the transition between the metastable and 
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FIG. 7: Dependence of Wif, on for various values of the anchoring energy parameter . It 
is shown that = w^^ ^ at wlp'^ = w^^ ^ and w,^ = at w^^^ = Wc (see Eq. (|59|) ). 

the equilibrium configurations involves penetrating the energy barrier A/ that separates the 
states with different half-turn numbers. This barrier can be seen from Figs. Efa) and Efa) 
where the free energ y fITlll is pl otted as a function of (3q and /3, respectively. 

Previous authors [ij, |l5l |l7| have supposed that the transitions may occur only if there 
is no energy barrier. Clearly, this assumption implies that the jumps take place at the end 
points of the stability intervals: /3o = = (1/2 + k)Ti ± w^, where the configuration with 
the half-turn number k becomes marginally stable {A^ = = 0) and loses its stability. 

These transitions are indicated by arrows in Figs. El El and El As is seen from Fig. El^c), in 
this case the upward and backward transitions: k ^ k + 1 and k + 1 k occur at different 
values of /3o: /?+ ^ and (3^~^^^ = (3^^ — 2w^, respectively. So, there are hysteresis loops in the 
response of CLC cell to the change in the free twisting number. 

We can now describe how the increase in the anchoring energy will affect the scheme of 
the transitions. To be specific, we consider the critical end point (3o = f3^\ so that for small 
anchoring energies with < tt/2 there are only two configurations: the marginally stable 
initial configuration with (3 = (3k and the equilibrium structure with (3 = (3k+i- In this case 
Eq. (|T!^ has at most two roots and the jumps will occur as transitions between the states 
which half-turn numbers differ by the unity, |AA;| = 1. 

At w,j) = tt/2, as shown in Fig. Eta), we have two marginally stable structures of equal 
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energy: and Pk+2- The newly formed structure Pk+2 being metastable at 7r/2 < < vr 
will have the free energy equal to the energy of the equilibrium configuration P^+i at = vr. 
So, as illustrated in Figs. Et^b) andEt^b), both transitions — >■ /c + 1 and k —>■ k + 2 are 
equiprobable and we have the bistability effect at the critical point under = t^. 

For vr < < 37r/2 there are three configurations: the initial configuration (3ki the 
metastable configuration [3k+i and the equilibrium structure /3fc+2- The configuration (3k+z 
being formed at = 37r/2 will define the equilibrium structure at 27i < < 3tt and so 
on. 

The general result for the critical point (3q = (3^^ can be summarized as follows. When 
(/ + 1/2)77 < < {1 + 3/2)71, in addition to the marginally stable configuration P^, there 
are / + 2 stable configurations with the half-turn numbers ranged from A; + 1 to A; + / + 2. 
The half-turn number of the equilibrium structure equals k + l + 1 under /vr < < (/ -|- l)7r. 

It immediately follows that the restriction < ti imposed by Belyakov and Kats [l^ on 
the anchoring strength requires the relaxation transitions to involve only two structures with 
I A A; I = 1. Our result shows that, when the anchoring parameter falls between /vr and 
(/ + l)7r, the half-turn number change is | AA;| = / + 1 for the transitions between marginally 
stable and the equilibrium states. Clearly, we can have the transitions with even AA; that 
involve topologically equivalent configurations with common parity j^, 23|. Such transitions 
may also be induced by the thermal director fluctuations without formation of defects even 
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if the anchoring is infinitely strong i^l- Though the mechanism under consideration is 



rather difi^erent, neglecting the director fluctuations can only be regarded as a zero-order 
approximation. 

Indeed, according to our remark at the end of Sec lIIIBj the expression for the fluctuation 
correlation function (jH^ implies its divergence upon reaching a marginally stable state where 
= 0. It means that taking the fluctuations into account will give the transition points 
located within the stability interval. This fluctuation induced shift may also suppress the 
hysteresis provided the mean square angle deviation ■>/ (0^), computed from Eq. at 
Pq = {1/2 + 771)71 and z = z' = d/2, and the anchoring energy parameter w^j, are of the same 
order. 



B. Asymmetric cells 

When the anchoring energies at the surfaces are different, wj^~^ 7^ 1^]+^ and e ^ 0, sm t>+ 
on the right hand side of Eq. (fT^ equals zero at P = 7r/2 + tt/c and, as demonstrated in 
Fig. Hfb), we have additional intersection points of the curves 7+(/3) and 7-(/5). It can be 
shown that the stability conditions are now deflned by Eq. and the twist parameters (3 
of the marginally stable conflgurations, where if^ = 0, can be computed as the stationary 
points of 7^. 

These points represent the local maxima and minima of 7^ and are located at /3 = 
{l/2 + k)7i± A/3. The equation for A/? is 

«;i+)(l + e)sin(A/?-earctan[6cotA/5]) = - '"f f f , (55) 

^ 1 + (e — Ij cos^ Ap 

where A/3 e [0,7r/2]. 

/^From the stability condition if^ > the values of (3 for stable conflgurations fall between 
the stationary points {k — l/2)7r + Af3 and {k + l/2)n — A/3, where the half-turn number k 
is the even (odd) integer depending on the parity. The function 7^ monotonically increases 
and varies from {k — l/2)7r — w,/) to {k + l/2)7r + Wfj) on the stability interval with the 
half-turn number k. The effective dimensionless anchoring parameter w^, as opposed to the 
case of equal anchoring energies with w^^^ = w^, is now given by 

Wff) = wlp~^ cos(A/3 — e arctan[e cot A/3]) — A/3. (56) 

Clearly, we can now follow the line of reasoning presented in Sec. IIV Al to flnd out the 
results concerning hysteresis loops and bistability effects that are quite similar to the case of 
equal anchoring strengths (see Figs.|Hfa)-(b)). There are, however, two important differences 
related to Eqs. and (jSHl)- 

If A/3 7^ 0, the intervals of /3 representing stable director conflgurations are separated by 
the gap of the length 2 A/3. The presence of this gap is illustrated in Figs. Efa) andEfb). 
Fig. ^b) shows the gap between stable branches of the dependence of the free energy on 
/3. The values of /3 within the gap represent unstable conflgurations and form the zone of 
"forbidden" states in the CLC cell. 

The graph of the A/3 vs w^^"* dependence is presented in Fig. IHl As expected, the gap is 

shown to disappear in the limit of equal energies, w^^'* = ^ Another and somewhat more 
interesting effect is that there is a small critical value of w[j^^ below which A/3 also vanishes. 
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In order to interpret this effect, we note that Eq. (fT^ with (3q = (1/2 + A;)7r has the only 
solution, f] = Pq, provided the azimuthal anchoring energy parameters meet the condition: 



2wi ^wi^^ 



\w\ ' -w\ 



Another form of this condition 

l4^^-4"^l>rf (58) 

implies that the difference between the azimuthal anchoring extrapolation lengths is larger 
than the cell thickness. For hybrid cells, similar inequality is known as the stability condition 
of homogeneous structures 2^ 2^, 1231] . 



In this case the gap disappears and the dependence of (3 on becomes continuous in the 
manner indicated in Figs. Efc) andlSfd). Given the value of "* the relation ()57j) yields the 
threshold value for the anchoring strength at the upper substrate: 

Wc = —j^. . (59) 

So, the jumps and the gap will vanish at w^^^ = Wc- Analogously, as illustrated in Fig. [71 

W(j, goes to zero at the critical point wlj^^ = Wc, while for large values of the dependence 

of on w^j^^ is approximately linear. 

In closing this section we discuss the limiting case of semi-strong anchoring where ^ — >■ 
oo. For this purpose we can combine the stabillity condition with Eq. ()18j) linking the 
free twisting parameter f3o and the twist parameter f3 of the helical structures chararcterized 
by the energy (fT^ . 

/^From the stabillity condition ()35|) the gap separating the stability intervals ranged be- 
tween {k - 1/2)tc + A(3 and {k + l/2)n - A/3 is given by 

\0, '< I- 

This result also follows from Eq. in the semi-strong limit e 1. Similarly, the expression 
for ()56|1 simplifies to the following form 

= wl^^hm{2A/3) - A(3 . (61) 

Eqs. dnni) and (fF)T|l explicitly show that the gap and the hysteresis loops both disappear 
len the c 
inequality 



when the anchoring strength is sufficiently small and 21^^^^ < 1. In other words, when the 



> d (62) 

is satisfied, the discontinuities turn out to be suppressed. 

Thus, in asymmetric CLC cells, the left hand side of Eqs. ()58p and ()62|) define the critical 
cell thickness below which the (3 vs curve becomes continuous. For semi-strong anchoring, 
this effect can be seen from the curves shown in Fig. |H1 
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V. DISCUSSION AND CONCLUSIONS 



In this paper we have studied how the pitch wavenumber q of the hehcal director con- 
figuration in the CLC cell depends on the free twisting number go at different anchoring 
conditions imposed by the cell substrates. It is found that this dependence is generally 
discontinuous and is characterized by the presence of hysteresis and bist ability. 

We have shown that asymmetry in the strengths of the director anchoring with the 
substrates introduces the following new effects: 

(a) the jump-like behaviour of the twist wavenumber is suppressed when the difference 
between the azimuthal anchoring extrapolation lengths is larger than the cell thickness; 

(b) the twist wavenumber intervals of locally stable configurations with adjacent numbers 
of helix half-turns are separated by gap in which the structures are unstable. 

Using stability analysis we have emphasized the idea that the instability mechanism 
behind the transitions between the helical structures is dominated by the in-plane director 
fiuctuations. These fiuctuations may render the structures unstable only if the anchoring 
energy is finite. 

In this case the height of the energy barriers separating the CLC states in the space of 
in-plane variables is determined by the surface potentials and is also finite. The mechanism 
can, therefore, be described as slippage of the director through the anchoring energy barrier. 
Interestingly, a similar mechanism can be expected be important in trying to extend the 
theory of Ref . jl^ , where shear- induced melting of smectic-A liquid crystals has been studied 
in the strong anchoring limit, to the case of weak anchoring. 

The part of our analysis presented in Sec. lIV Al relies on the assumption that the transition 
between configurations with different half-turn numbers occurs when the initial structure 
loses its stability, so that pitch wavenumber is no longer a local minimum of the free energy 
surface. The result is that the stronger the anchoring, the larger the change of the half-turn 
number (and of the twist wavenumber) needed to reach the equilibrium state. So, whichever 
mechanism of relaxation is assumed, the metastable states certainly play an important part 
in the problem when the anchoring is not too weak. 

It was recently shown by Bisi et. al jl^l that the instability mechanism in twisted nematics 
may involve the so-called eigenvalue exchange configurations (isL |29| . These configurations 
and the tilted structures are, however, of minor importance for the director slippage induced 
instability. They may be important outside the parameter regime considered here, and we 
will discuss alternative mechanisms in more detail elsewhere. 

The djTiamics of the transitions is well beyond the scope of this paper. Despite some very 
recent progress it still remains a challenge to develop a tractable theory that properly 
account for director fiuctuations, hydrodynamic modes and defect formation. Simultane- 
ously we have seen at the end of Sec. IIV Al that fiuctuation effects can be estimated by using 
the expression for the correlation functions given in Sec. IIII CI But in order to take the 
fluctuations into consideration a systematic treatment is required. 
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APPENDIX A: FLUCTUATION SPECTRUM AND STABILITY CONDITIONS 

In this appendix we comment on the eigenvalue problem written in a form similar to 

Eqs. dSHj-dsni) 

[dl + A] Xx{t) = 0, (Al) 
±drXx + w±X^ = 0, (A2) 

- T=±l 

where A is the eigenvalue and X\{t) is the eigenf unction. Our task is to derive the conditions 
which ensure positive definiteness of the eigenvalues. 

To this end we consider the case of negative eigenvalues with A = — and substitute the 
general solution of Eq. ()A1|) 

Xx{t) = Ax sinh kt + Bx cosh kt (A3) 



into the boundary conditions ()A2|) . This yields a homogeneous system of two linear algebraic 
equations for the intergration constants Ax and Bx- The system can be written in matrix 
form as follows 

"■{I 

where 



_ I K cosh K + w+ sinh k k sinh k + Wj^ cosh k , . 
' cosh K — w_ sinh K k sinh k + cosh k ' ' 



Non-zero solutions of Eq. ()A4|) exist only if the determinant of the coefficient matrix H 
vanishes, detH = 0. For the matrix ()A5|) . this yields a transcendental equation 



sinh 2k + {w+ + wJ)Kcosh. 2k + Wj^w_ sinh 2k = (A6) 

whose roots determine the negative eigenvalues through the relation A = — k^. 
Eq. flA6|l can be conveniently recast into the form 

/(x) = + ax coth X = b, (A7) 

where x = 2k, a = 2(u;_|_ + w_) and b = —Aw^w^. It is now not difficult to see that the 
inequahty 

b < min f(x) = minfa;^ + ax cothx) (A8) 

x>0 x>0 

provides the condition for the eigenvalues to be positive. 

Eq. flA8|) combined with the relation b > — can be analyzed using elementary meth- 
ods. The final result 

b< a, a > (A9) 
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immediately leads to the inequalities 

w+ + w_ + 2w+w_ > 0, w+ + w_ > 0, (AlO) 

which ensure positive definiteness of the spectrum. In the semi-strong limit, w_ — > oo, the 
conditions (jAlOjl assume the following simplified form: 

l + 2w+> 0. (All) 

Our final remark is that changing k to m in Eq. ()A6jl gives the equation 

{w-f^W- — K^) sin 2k + {w^ + wJ)kcos 2k = 0. (A12) 

The roots of this equation determine positive eigenvalues A = k^. In the strong anchoring 
limit, w± —>■ oo, Eq. ()A12|) takes the well known form sin 2k = leading to the relation 
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